Magnetic Properties of the S=| Quasi-One-Dimensional Antiferromagnet CaCu203 
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We report single crystal growth and magnetic susceptibility and neutron diflraction studies of the 
S=i quasi-one-dimensional antiferromagnet CaCu2 03. The structure of this material is similar to 
that of the prototype two-leg spin-ladder compound SrCu2 03. However, the Cu-O-Cu bond angle in 
the ladder rungs in CaCu2 03 is equal to 123°, and therefore the magnetic interaction along the rung 
is expected to be much weaker in this material. At high temperatures, the magnetic susceptibility 
of CaCu2 03 can be decomposed into a contribution from one-dimensional antiferromagnetic chains 
or finite size chain segments together with a weak Curie contribution. The intrachain magnetic 
exchange constant ,J\\ , determined from the magnetic susceptibility measurements, is 2000±300 K. 
CaCu203 undergoes a Neel transition at Tjv=25 K with ordering wavevector of (0.429(5), |, |). 
The magnetic structure is incommensurate in the direction of the frustrated interchain interaction. 
Weak commensurate magnetic Bragg peaks with the reduced wavevector (i, i, i) are also ob- 
served below Tjv. Application of a magnetic field induces a metamagnetic transition at which the 
incommensurability of the magnetic structure is substantially reduced. The material possesses only 
short-range magnetic order above the transition field. 

PACS numbers: 75.25.-hz, 75.40. Cx, 75.50.Ee 



I. INTRODUCTION 

Low-dimensional quantum spin systems exhibit a va- 
riety of intriguing properties originating from the low 
dimensionality and quantum fluctuations. Quantum ef- 
fects are most pronounced in S=^ systems. The one- 
dimensional (ID) S=i nearest neighbor Heisenberg anti- 
ferromagnetic (AF) chain and its two-dimensional (2D) 
analog, the square lattice antiferromagnet, have been ex- 
tensively investigated both theoretically and experimen- 
tally. The properties of these systems turn out to be very 
different: the nearest-neighbor Heisenberg AF chain ex- 
hibits a critical ground state ||l|, while true long range 
order takes place at temperature T=0 in the 2D sys- 
tem Pj. These observations have motivated extensive 
research activity that has resulted in the discovery of 
magnetic systems consisting of a finite number of inter- 
acting ID chains which interpolate between the ID and 
2D cases. 

The most notable examples of such systems are pla- 
nar arrays of strongly coupled n chains, the so-called 
n-legged spin ladders |^. 8=^ n-legged spin ladders ex- 
hibit fundamentally different properties for even and odd 



n. Even-legged ladders exhibit a spin liquid ground state 
with a gap in the magnetic excitation spectrum, while 
odd-legged ladders exhibit a ground state isomorphous 
to that of the ID chain. Several good physical realiza- 
tions of these systems, notably Srin-iCum+i02m, and 
Sri4Cu2404i, have been found, and the main predictions 
for the nature of the ground state have been confirmed 
experimentally 0]. However, our understanding of these 
intriguing materials is still incomplete from both the the- 
oretical and experimental points of view. 

Adding to the complexity of the picture, frustrated 
magnetic interactions and disorder are often present in 
real-life materials and result in new and unexpected phe- 
nomena. In the cuprate spin-ladder systems mentioned 
above, the interladder Heisenberg interactions are indeed 
essentially frustrated. That is, classically they cancel ex- 
actly, and therefore weak magnetic interactions resulting 
from magnetic anisotropy or from minute lattice distor- 
tions should play an important role in these compounds. 
The effects of frustration and disorder together with the 
effects of anisotropic interactions on the properties of low- 
dimensional materials provide a rich field for both the- 
oretical and experimental work. A substantial amount 
of research activity is currently being devoted to these 



subjects. 

In this paper, we report single crystal growth and mag- 
netic susceptibility and neutron diffraction studies of the 
S=i quasi-lD magnet CaCu203. The structure of this 
material, shown in Fig. m, is similar to that of the pro- 
totype two- leg spin-ladder compound SrCu203. It con- 
sists of an array of ladder-like structures with quasi- ID 
copper-oxide chains extending along the crystallographic 
b direction. However, while the Cu-O-Cu bond angle in 
the rungs is 180° in the Sr compound M, it equals 123° in 
CaCu203; colloquially, the ladder is "buckled". There- 
fore, the interchain interaction in the buckled ladders in 
CaCu203 should be substantially weaker than that in 
SrCu203. At high temperatures, the temperature de- 
pendent part of the magnetic susceptibility of CaCu203 
can be described as the sum of a contribution from one- 
dimensional antiferromagnetic chains or finite size chain 
segments, and a weak Curie contribution. The intra- 
chain magnetic exchange constant Ju , determined from 
the slope of the high-temperature susceptibility curve, is 
J|[=2000±300 K. Similar values have been previously re- 
ported for the related quasi-lD compounds Sr2Cu03 and 
SrCu02 i- 0. 

CaCu203 undergoes a Neel transition to 3D long-range 
order at Tjv ~25 K. The magnetic ordering occurs with 
a three-dimensional-like critical exponent (3 = 0.3 — 0.4. 
The ordering wavevector is (0.429(5), i, i). The mag- 
netic structure is incommensurate in the direction of the 
frustrated interchain interaction. Weak commensurate 
magnetic Bragg peaks with the reduced wavevector (i, 
h, h) are also observed below Tjy. The experimental data 
are consistent with a spiral (or possibly conical) magnetic 
structure with spins rotating in the ac crystallographic 
plane. The low-temperature magnitude of the ordered 
magnetic moment determined using this model magnetic 
structure is 0.2±0.07 fiB- To explain the observed mag- 
netic structure, either the presence of anisotropic mag- 
netic interactions or the existence of a lattice distortion 
lifting the frustration of the intcrladder Heisenberg inter- 
action seems to be required. The application of a mag- 
netic field induces a metamagnetic transition in which the 
incommensurability of the magnetic structure is substan- 
tially reduced. The material possesses only short-range 
magnetic order above the transition field. 

The format of this paper is as follows. In Section II, we 
describe the single crystal growth and the sample char- 
acterization, together with experimental details of the 
magnetic susceptibility and neutron diffraction measure- 
ments. Section III presents the results and analysis of the 
magnetic susceptibility measurements. In Section IV, we 
discuss neutron diffraction measurements of the crystal 
structure, zero-field magnetic ordering, and also the ef- 
fects of a magnetic field. Finally, in Section V we give a 
summary and conclusions. 



II. EXPERIMENTAL DETAILS 

CaCu203 has an orthorhombic lattice, space group 
Pmmn (Ref. @), with a=9.949 A, 6=4.078 A and 
c=3.460 A at T=10 K. Single crystals of CaCu203 
have been grown using the traveling solvent floating 
zone (TSFZ) method with CuO as a flux. To prepare 
CaCu203 powder, a mixture of CaCOs and CuO pow- 
ders was calcined at 850 C for 12 hours with subsequent 
heating at 995 C in air for 3 days with frequent grind- 
ings. Only small amounts of Ca2Cu03 and CuO impuri- 
ties have been found in the as-prepared powder. Oxygen 
or inert gas annealing at 995 C does not result in any 
signiflcant magnetic or structural change, which implies 
that the oxygen content is close to being stoichiometric. 

The growth is carried out in a four-lamp optical FZ 
furnace built by Crystal Systems Inc. The temperature 
range, in which crystals of the required composition are 
stable, 975-1020°C, is extremely narrow Q, requiring 
careful temperature control. We use as-prepared powder 
to make feed rods. To reduce the radial thermal gradi- 
ent of the hot zone, feed rods with diameter less than 5 
mm are used. By careful observation of the amount of 
extra feed rod that dissolves into the molten zone, we 
estimate that 17.8% CaO / 82.2% CuO is the optimal 
initial concentration of the flux rod. Stable growth is 
achieved using a pulling rate of 0.8 mm/hour, a 35 rpm 
rotation rate, and an ambient gas flow of oxygen (220 
ml/min) and argon (20 ml/min). The growth direction 
in the obtained crystals is approximately parallel to the 
copper-oxide chain (010) direction. 

Magnetic susceptibility measurements in the tempera- 
ture interval 2-800 K were carried out using a commercial 
Quantum Design SQUID magnetometer equipped with 
a furnace attachment. To perform measurements with 
the magnetic field oriented along chosen crystallographic 
axes, the crystals were aligned using x-ray diffraction 
techniques. 

For single crystal neutron diffraction experiments, two 
single grain crystals of CaCu203 with mosaicities of 
r-^ 0.2° from two different batches were chosen. We refer 
to them as sample 1 and sample 2. The samples were of 
cylindrical shape, ^--^4 mm in diameter and ~20 mm in 
height. Neutron powder diffraction measurements per- 
formed on powder samples prepared from various other 
single crystals showed that up to 7% by weight of the 
impurity phase Ca2Cu03 may be present. However, 
no single-crystalline phase other than CaCu203 was de- 
tected in the investigated samples. Therefore, any impu- 
rity phase, if present, was in the form of scattered powder 
inclusions. 

The neutron diffraction measurements were carried out 
at the reactor of the National Institute of Standards 
and Technology (NIST). The powder diffraction experi- 
ment was carried out at the BT-1 high- resolution powder 



diffractometer. The data were collected in the 20-range 
3° — 165° using a neutron beam of wavelength 1.5401 A 
produced by a copper (311) monochromator. The coUi- 
mations used were 15'-20'-S-7'. Data were collected at 
295, 40, and 10 K. Crystal structure refinements were 
carried out using the GSAS program [Q. Single crys- 
tal measurements were performed at the BT-2 and BT-9 
triple axis spectrometers, using neutrons of energy 14.7 
and 30.5 meV for the measurements with and without 
applied magnetic field, respectively. Pyrolytic graphite 
(002) monochromator and analyzer crystals were used. 
To reduce higher energy harmonics present in the beam, 
pyrolytic graphite filters were placed before and after the 
sample. The collimations used were 60'-60'-S-80'-80' and 
60'-40'-S-40'-40' for the measurements taken in zero and 
applied magnetic field, respectively. Zero-field experi- 
ments were performed using either a closed cycle helium 
refrigerator with base temperature lOK, or a pumped 
helium cryostat which allowed one to perform measure- 
ments down to T=1.5K. Data were taken in the (h, k, k), 
(h, k, 3k), (h, 3k, k), and (h, 5k, k) scattering planes. A 
9 Tesla vertical field superconducting magnet was used 
for the measurements in an applied field, and the experi- 
ment in a field was performed in the (h, 3k, k) scattering 
plane. 



III. TEMPERATURE-DEPENDENT 
SUSCEPTIBILITY 

Figure shows the temperature dependence of the 
magnetic susceptibility of CaCu203 for T<300 K mea- 
sured in a magnetic field of 1 kG aligned with 
the three major crystallographic axes. Temperature- 
independent contributions of 2.23-10"'* cm'^/mol Cu, 
1.73-10"'* cmVmol Cu, and 2.37-10"'' cm^/mol Cu were 
subtracted from the data of Fig. || for the magnetic field 
aligned with the a, b, and c crystallographic axes, re- 
spectively. The sample undergoes a Neel-like transition 
at Tat ~25 K. The anisotropy of the susceptibility be- 
low Tn indicates that spins are mostly confined to the 
ac crystallographic plane at low temperatures. For tem- 
peratures 50 K < T < 300K, the temperature depen- 
dence of the magnetic susceptibility is well described by 
Curie 1/T behavior (see the inset in Fig. g). The effec- 
tive concentration of Cu^^ free spins extracted from the 
corresponding Curie constant amounts to 3% of all the 
Cu atoms present in the sample. It is impossible, how- 
ever, to ascribe the Curie term to extraneous impurities 
since the Curie contribution is dramatically decreased be- 
low the Neel temperature, and therefore the "free" spins 
responsible for the Curie term participate in the low- 
temperature magnetic ordering. Specifically, the data of 
Fig. imply that the effective concentration of Cu'^^ 
free spins present in the sample below T=25 K is less 
than 1%. As shown in Section IV, our neutron studies 



show that at T^r CaCu203 undergoes a transition to a 
magnetically ordered state with a relatively large value 
of the ordered magnetic moment, for a quasi-lD S = 1/2 
antiferromagnetic system. This, in turn, means that vir- 
tually all of the Cu^+ spins participate in the magnetic 
ordering. 

From Fig. one sees that for TkT^t the material 
undergoes a field induced transition when H is parallel 
to the crystallographic a axis. The transition results in 
a small increase of the magnetic moment at H~ 3T at 
low temperature, and exhibits hysteresis on ramping the 
magnetic field up and down for T< 20 K. For H parallel 
to the c axis there is no jump in the moment but the slope 
of M(H) increases at approximately the same field. (See 
Fig. y|b.) The inset in Fig. ^ shows the temperature de- 
pendence of the transition field He for the magnetic field 
parallel to the a crystallographic axis. To determine He 
we have averaged the maxima in the field derivative of 
the magnetization for increasing and decreasing the mag- 
netic field. The error bars in the inset reflect the width of 
the hysteresis region for T<20 K; at higher temperatures, 
error bars correspond to the with of the peak in dM/dH. 
The transition shown in Fig. exhibits some features 
characteristic of a spin-fiop transition in a magnetic sys- 
tem in which spins are confined to the ac crystallographic 
plane at zero field, but are not parallel to either a or c. In 
particular, the transition does not occur when the mag- 
netic field is parallel to the b axis. The transition field He 
varies with temperature in the way characteristic to the 
spin-flop fleld of a typical anisotropic antiferromagnet [|| . 
For H parallel to c the transition appears to be contin- 
uous even at the lowest temperatures. Such continuous 
transitions are known to occur in coUinear antiferromag- 
nets when the easy axis is not parallel to the applied field 
Q, and also in compounds with non-coUinear magnetic 
structures. In Section IV we show that the magnetic 
structure of CaCu203 is, in fact, non-coUinear. 

The dominant magnetic interaction in CaCu203 is the 
intrachain coupling Ju which, by analogy to Sr2Cu03 and 
SrCu02, is expected to be of the order of 2000 K. How- 
ever, the spin gap is proportional to the leg coupling J±. 
In SrCu203 J± ~ ^J\\ jl^l and the spin gap is equal to 
420 K ||l^. Because of the reduced value of J± in the 
buckled ladders of CaCu203, we expect that the mate- 
rial will have a smaller gap and will exhibit ID behavior 
characteristic of uncoupled spin-chains above T~300-400 
K [l4| . At lower temperatures, three dimensional mag- 
netic interactions are, in general, not negligible, and the 
magnetic susceptibility can only be calculated accurately 
if these interactions are taken into account. To estimate 
the numerical value of the intrachain coupling Jii , we have 
used two approaches which are described below. 

Figure H shows high temperature susceptibility mea- 
sured in a field of 1 Tesla approximately aligned with the 
spin-chain direction. The Curie term is subtracted from 
these data. The sample starts to decompose at T~700 K, 



and, as a result, the susceptibility abruptly decreases for 
T>700 K. To analyze the contribution to the magnetic 
susceptibility from the quasi-lD S=^ system Xspin, we 
write X ^s 



X = XCn 



Xc 



Xv 



" Xspin 



(1) 



where Xcore is the diamagnetic contribution from the 
ionic cores, and Xw is the Van Vleck paramagnetic term. 
Xcore and Xvv are independent of temperature. Equation 
is phenomenological because, as we have mentioned 
above, the low-temperature Curie term cannot be as- 
cribed to noninteracting extraneous impurities, and thus, 
at this stage, Eq. Q lacks a physical explanation. The 
high-temperature data of Fig. ^ were fitted to the theo- 
retical results for x(T) for a ID 8=^ chain obtained us- 
ing the Bethe ansatz method ||l^. The solid and dashed 
curves in Fig. ^ are Bethe ansatz results for three dif- 
ferent values of Jii with a temperature independent term 
representing Xcore + Xf« added. Evidently, above T~300- 
400 K and below the sample decomposition temperature, 
Xspin (T) is well described by the theoretical curve for a 
ID S=i chain with J=1950±300 K. This value is in good 
agreement with the results for the related quasi- ID chain 
compounds SraCuOg, SrCuOz, and SrCuaOg [|-0. Us- 
ing the value Xcore=-3.3-10~^ cmVmol Cu2+ (Ref. [|§), 
we obtain Xw — ^-O-W^'^ cm'^/mol Cu. This result for Xw 
is somewhat larger than the previously reported values 
of x™=2-8-10-5 cmVmol for SrzCuOs and SrCuOz @. 
We note that these results are virtually independent of 
the details of the low-temperature fits that have been 
used to extract the XCurie term. 

To explain the observed temperature dependent mag- 
netic susceptibility at a microscopic level, we consider 
a model in which impurities or structural imperfections 
break the quasi- ID magnetic chains into weakly inter- 
acting segments of finite size. We have calculated the 
susceptibility of finite spin-chain segments using quan- 
tum Monte Carlo simulations utilizing the loop cluster 
algorithm ||l^,|l^. The only parameters in these calcu- 
lations are the size of the chain segments, and the in- 
trachain magnetic coupling J||. We use 3=2.1 for the 
Lande factor, which is a typical value for Cu^"*" in insula- 
tor cuprate compounds [Q. Chain segments with even 
and odd number of spins exhibit qualitatively different 
behavior. In our calculations, we assume that identi- 
cal numbers of "even" and "odd" segments are present. 
We find that this model reproduces the Curie-like behav- 
ior of the magnetic susceptibility at low temperatures, 
as well as the gradual increase of the susceptibility with 
temperature at high temperatures. However, to achieve 
satisfactory agreement with the experimental data, it is 
necessary to add a small additional Curie term to the 
calculated susceptibility. The magnetic susceptibility in 
this model is, therefore, given by the same equation as 
the phenomenological model discussed above (Eq. |^) in 



which Xspin now stands for the magnetic susceptibility of 
the finite ID chain segments. 

In Fig. ^ we show the CaCu203 magnetic suscep- 
tibility data together with the results of Monte Carlo 
calculations for the susceptibility of equal numbers of 
spin-chains consisting of segments with 13 or 14 spins 
in every segment. We use J|| = 1940 K for the intrachain 
exchange. A Curie term corresponding to a 1% effec- 
tive concentration of Cu'^^ free spins, as well as Van 
Vleck x-u« =9.4-10"^ cm-^/mol Cu and ionic core Xcore=- 
3.3-10"^ cm'^/mol Cu contributions are added. Calcu- 
lations for the finite-segment model with segments con- 
taining more than 20 spins, or less than 10 spins, are 
incompatible with the experimental data. Thus, in the 
context of this model, a substantial number of weak links 
in the magnetic chains is required to explain the ob- 
served magnetic susceptibility. Such weak links can occur 
due to presence of oxygen vacancies or extrinsic impuri- 
ties. Whether or not the actual concentration of defects 
in the crystal structure of our samples can explain the 
large number of weak links required by the finite-segment 
model will be the subject of future work. 

The results of our calculations exhibit systematic devi- 
ations from the experimental data for temperatures be- 
low T=300 K. This indicates that, as expected, inter- 
chain interactions are not negligible at these tempera- 
tures. However, our Monte Carlo calculations, that were 
carried out in absolute units with a small number of ad- 
justable parameters for a simple finite-segment model, 
are clearly in good qualitative agreement with the ex- 
perimental data. Therefore, we believe that the finite- 
segment model provides a plausible explanation for the 
unusual temperature dependence of the magnetic suscep- 
tibility, including the presence of the Curie-like term at 
temperatures higher than T^v- Since all of the copper 
spins are expected to participate in the low-temperature 
3D magnetic ordering, this term should disappear below 
the Neel temperature, in agreement with the experimen- 
tal results. 

At temperatures higher than T=300 K, the agreement 
of the finite-segment model with the experimental data 
is satisfactory (see the inset in Fig. |5|) . The value of the 
intrachain coupling J|| = 2000±500 K extracted from the 
high-temperature fits to this model is insensitive to both 
the size of the chain segments and to the magnitude of the 
added Curie term. Therefore, the two models described 
above give the same value of Jh within the errors. 



IV. NEUTRON DIFFRACTION 

A. Zero-Field Measurements 

The structure of CaCu203 is shown in Fig. n^. The cor- 
responding atomic coordinates and temperature param- 
eters determined in the neutron powder diffraction ex- 



periment at T=10 K (Rp=3.53%, Ru,p=4.45%) are listed 
in Table |. The structure contains copper oxide chains 
running along the crystallographic b direction; the chains 
form ladder-like pairs as shown in Fig. n^. As we have 
mentioned above, the main structural difference between 
CaCu203 and its Sr analog is that in CaCu203 the Cu- 
0-Cu angle in the rungs of the ladders is 123°, while 
the corresponding angle equals 180° in SrCu203. As 
in SrCu203, the ladders are coupled to each other via 
~90° Cu-O-Cu bonds in the a direction. In the crystallo- 
graphic c direction, the copper-oxide ladders are stacked 
on top of each other with Ca atoms in between. 

CaCu203 undergoes a 3D Neel transition at T^ ~25 
K. Figure |6| shows an elastic neutron diffraction scan 
along the (h, 0.5, 0.5) direction at T=12 K. (The 
wavevector transfers are quoted in reciprocal lattice 
units, r.l.u.) Below T=25 K, superlattice Bragg peaks 
with reduced wavevector (0.429(5), 0.5, 0.5) are detected. 
Note that, within the accuracy of our experiment, these 
peaks can be described as commensurate peaks at the 
(3/7, 0.5, 0.5) position. In this case, the magnetic unit 
cell is seven times larger than the chemical unit cell in 
the crystallographic a direction. Much weaker scatter- 
ing is also found at the (0.5, 0.5, 0.5) positions. The 
latter peaks are not due to nuclear multiple scattering 
effects since they are temperature dependent and disap- 
pear above Tjv (see Fig. uh. The intensity ratio of the 
(0.5, 0.5, 0.5) peak to the (0.429, 0.5, 0.5) peak is sam- 
ple dependent, and specifically, in sample 2 this ratio is 
approximately two times smaller than that in sample 1. 

The data of Fig. o and the corresponding data taken at 
various temperatures have been fitted to three resolution- 
limited Gaussian peaks. The position of the central peak 
was fixed at the commensurate value of (1.5, 0.5, 0.5). 
The temperature dependence of the magnetic Bragg peak 
intensities is shown in Fig. 0. The critical exponent (3 
and the Neel temperature Tjv are extracted by fitting the 
data to the power law / ~ (1 — T/Tn)^^ ■ For samples 1 
and 2 the parameters are found to be Tjv=25.7(0.2) and 
25.4(0.2) and /3=0.35(0.03) and 0.44(0.04), respectively. 
The large values of f3 indicate that the transition has a 
three dimensional Heisenberg character. The differences 
between the samples are likely due to slight variations 
in their composition. Despite these differences, the posi- 
tions of the incommensurate peaks are the same in the 
both samples. The positions of these peaks are temper- 
ature independent. 

Neutron diffraction data were also collected in the (h, 
k, 3k), (h, 3k, k), and (h, 5k, k) reciprocal space zones. 
Altogether, the integrated intensities for 27 magnetic 
Bragg reflections were measured at T=12 K. These data 
were used for the determination of the low-temperature 
magnetic structure described in the next section. 



B. Zero-Field Magnetic Structure 

In this section we present a simple model magnetic 
structure consistent with our experimental observations. 
To construct a plausible model, we first consider the mag- 
netic interactions known to be present in CaCu203. The 
magnetic interaction is strongest along the copper-oxide 
chains. As shown in Section III, the intrachain Heisen- 
berg exchange constant is J||=2000±300 K. According to 
the Goodenough-Kanamori- Anderson rules [|2[ , the rung 
exchange constant J± is greatly reduced compared to J\\ 
because of the reduced Cu-O-Cu angle in the rungs. Due 
to the diminished rung coupling, the energy gained upon 
the transition to the gapped spin-liquid state is also di- 
minished since this gain is proportional to the gap value, 
which goes to zero in the limit of decoupled chains [|j. 
Therefore, the three-dimensional magnetic interactions 
in CaCu203 are relatively more important than those in 
SrCu203. The fact that the former compound undergoes 
a 3D Neel transition at TAr=25 K while the latter stays in 
the non-magnetic state down to the lowest temperatures 
||l3| is consistent with this suggestion. 

However, since the Cu-O-Cu angle is different 
from 180°, in addition to the Heisenberg coupling a 
Dzyaloshinsky-Moriya (DM) interaction is also present 
in the rungs of the ladder. This interaction is of the form 
||l9| Hdm = D • Sj X S-, . The magnitude of the DM vector 
D can be estimated [Q as D^ —^Jj_^, where A^ is the 
shift of the gyromagnetic ratio g from the g-value for a 
free electron, and $ is the deviation of the Cu-O-Cu rung 
angle from 180°. Since Ag/g is typically of the order of 
0.1, and <I> ~1, D may be as large as several meV. The 
direction of the DM vector is defined by the lattice sym- 
metry |2^. In CaCu203, it points along the 6-axis; the 
DM vectors for the two ladders belonging to the same 
chemical unit cell (see Fig. |l|b,c) are antiparallel. The 
DM interaction is anisotropic with the anisotropy energy 
of the order of D^/2J_l. In the present case, it favors 
configurations with spins confined to the ac-plane. 

The strongest interladder interaction is likely to be 
along the c axis, in the stacking direction of the lad- 
ders. As pointed out by Greven and Birgeneau |2^, the 
exchange constant along this direction, Jc, should be of 
the order of 10 meV in SrCu203, and because the geom- 
etry is similar, it is likely of the same order of magnitude 
in CaCu203. The interladder interaction in the a direc- 
tion is more complex. The Heisenberg interaction in this 
direction is ferromagnetic and relatively weak (J^ ~10 
meV) ll^. Moreover, because of the extremely strong 
AF intrachain interaction Jy , the magnetic coupling in 
the a direction is essentially frustrated. Specifically, in 
the classical limit the a-direction interladder interaction 
cancels exactly. Therefore, to a first approximation, the 
magnetic system consists of an array of 2D double planes 
in which spins are coupled by Heisenberg interactions 



with antiferromagnetic coupling constants J||, Jc^ and 
Jj_. In the normal direction (crystallographic a axis), 
the interplane Heisenberg interaction is frustrated, and 
the planes are only weakly coupled. A non-zero interpla- 
nar magnetic coupling, which is required for the 3D mag- 
netic ordering to occur, can result from weak anisotropic 
interactions 123^24], such as the pseudo-dipolar interac- 
tion. Another possibility is that an appropriate lattice 
deformation, small enough to be undetected in our pow- 
der measurements, lifts the exact cancellation of the in- 
terladder Heisenberg interaction, thereby providing the 
necessary interplanar coupling |^,^ . 

In the following refinement, we address the magnetic 
structure that gives rise to the strong incommensurate 
reflections; the origin of the weak commensurate peaks 
will be briefly discussed at the end of this section. Mag- 
netic susceptibility measurements discussed in Section 
III indicate that spins are mostly conflned to the ac- 
crystallographic plane, consistent with the geometry of 
the DM interaction. Therefore, we start with trial struc- 
tures obeying this restriction and then remove it. Since 
the unit cell is doubled in the b and c directions, the mag- 
netic structure is of a simple antiferromagnetic type along 
these axes. At the moment, it is unclear how the inter- 
ladder interactions give rise to the observed incommen- 
surability along the a axis; a detailed theoretical model 
is needed to answer this question. In the absence of a 
theory, we have considered two simple phenomenological 
models for the incommensurate spin structure: a spin 
density wave model, and a spiral model. We have found 
that the former does not reproduce the experimentally 
observed magnetic peak intensities. As discussed below, 
the spiral model is, on the other hand, in reasonable 
agreement with the experimental data. 

In the spiral model (Fig. |l|c), we assume that the 
spins rotate by a constant angle x from ladder to lad- 
der. The net interladder interaction along a is expected 
to be much weaker than the interaction across the rungs, 
and therefore the relative orientation of the rung spins 
is defined by the latter interactions. Since both Heisen- 
berg and DM interactions are present in the rungs, the 
angle between the rung spins should deviate from 180°. 
We denote this angle as 77. In the classical approxima- 
tion, tan 77 = D/ J±_ ^ —^^. The sign of 77 changes from 
ladder to ladder together with the direction of the DM 
vector (see Fig. ||c). 

At low temperatures, the integrated intensity of the 
magnetic Bragg peaks is given by [M 






sin(26l) 



Q/3 



fc„fc^) 5^(5,'?) (5f)exp{zk. (1-1')} 
w 

(2) 



indices x,y and z, and F(k) is the Cu^+ magnetic form 
factor. The sum is taken over all copper atoms in the 
magnetic unit cell. Since the reduced wavevector of the 
magnetically ordered state equals, within the accuracy of 
our experiment, (3/7, 0.5, 0.5), we have assumed that the 
actual magnetic unit cell is 7a x 26 x 2c, expressed in terms 
of the chemical unit cell parameters. The anisotropic 
3(^3:2 _j,2 magnetic form factor was calculated using for- 
mulae given in Ref. |^ . The period of the lattice mod- 
ulation fixes X ^t 37r/7-|-7rn, where n is an integer. 

We use our spiral model spin structure and Eq. || to 
fit the experimental data which consists of 27 magnetic 
Bragg reflections collected in the (h, k, k), (h, k, 3k), (h, 
3k, k), and (h, 5k, k) zones. In addition to 77 and the spin 
magnitude S, the tilt of the plane of spin rotation with 
respect to the crystallographic axes is varied. The best 
results are achieved for the spins rotating in the ac-plane 
with 77 = 160° ± 20°. The crystallographic reliability in- 
dex, deflned as R=X] Wl observed- Vl calAI Z) ^observed, 

is R=0.18 (see Table |l|). The magnitude of the ordered 
magnetic moment calculated in this model by comparing 
the magnetic peak intensities with the intensity of the 
weak structural Bragg peak (5,1,1) is 0.2 ± .07/i_b. 

Despite the relatively large value of R, the simple 
model described above reproduces the qualitative behav- 
ior of the experimental data, such as the high intensity 
of the peaks with h=1.43 and 4.43 as compared to those 
with h=0.43, 2.43, and 3.43, together with the growth of 
the peak intensity with the projection of the scattering 
unit-vector on the 6 axis. The value of 77 is consistent 
with the rough estimate of the DM constant D given 
above, and the fitted magnitude of the ordered magnetic 
moment is typical for quasi-lD S=2 niagnetic systems. 
Taking into account the simplicity of the model, these 
results are quite satisfactory. However, our model is in 
all likelihood oversimplified. It does not, for example, ex- 
plain the anisotropy of the magnetic susceptibility in the 
ac plane (see Fig. Q). We have considered elliptic mag- 
netic structures that could account for this anisotropy. 
Such structures, however, produce fits of similar quality 
to those obtained using the circular spiral model. Clearly, 
more theoretical and experimental work is required to 
elucidate the magnetic structure of CaCu203 in detail. 

The weak commensurate scattering at the (0.5, 0.5, 
0.5) position can have at least two different origins. Cant- 
ing of the spins out of the ac-plane with canting angle 
of the order of 10° and with the appropriate periodicity 
is one possible scenario. In this case, a conical mag- 
netic structure is realized. Alternatively, commensurate 
and incommensurate components can originate from spa- 
tially distinct parts of the sample, possibly differentiated 
by slight compositional variations. 



where 20 is the scattering angle, k is the scattering vec- 
tor, 1 denotes atomic positions, a, (3 stand for the space 



C. Effects of a Magnetic Field 

The application of a magnetic field has a dramatic 
effect on the low-temperature magnetic structure in 
CaCu203. According to the magnetic susceptibility mea- 
surements of Section III, application of a magnetic field 
along the a or c axis results in a transition at a mag- 
netic field of approximately 3.5 Tesla at low tempera- 
tures. Taking into account the zero-field magnetic struc- 
ture discussed in the previous section, it is natural to 
expect that, in agreement with the magnetization data 
of Section Iff, application of a magnetic field along the a 
or c axis may result in a spin-flop-like transition. In gen- 
eral, the periodicity of the magnetic structure is expected 
to change in this transition. 

To investigate the effects of a magnetic field, we have 
performed elastic neutron diffraction scans in the vicinity 
of the (0.5, 1.5, 0.5) reciprocal lattice position in sample 
2. The horizontal scattering plane coincided with the 
(h, 3k, k) zone, and the field was perpendicular to the 
scattering plane. In this geometry, the magnetic field has 
no a-axis component, and the angle between the field and 
the c-axis is 22°. 

The diffraction scans taken at T=10 K in various mag- 
netic fields are shown in Fig. |^. As may be seen in Fig. 
|8|, application of a magnetic field results in a substan- 
tial reduction of the incommensurability of the magnetic 
structure. However, even in a field of 8 Tesla, the mag- 
netic structure in the a-direction is still incommensurate. 
An unexpected result is the measurable broadening of the 
diffraction peaks with the application of a magnetic field. 
This broadening indicates that the long-range magnetic 
order characteristic of the zero-field magnetic structure 
is destroyed at high magnetic fields. 

The data of Fig. have been fitted to the sum of 
three Gaussian peaks. The results of the fits are shown 
as solid lines in Fig. |8[ Fig. ^ shows the magnetic field 
dependence of the intensities of all the three peaks, the 
separation between the incommensurate peaks, and the 
full peak width at T=10 K. In agreement with the sus- 
ceptibility data for the magnetic field parallel to the c 
axis, the transition begins at H~3 Tesla and is not com- 
plete at H=5 Tesla. In fact, it is evidently not complete 
even at H=8 Tesla, as the data of Fig. g show. The 
redistribution of the scattered intensity between the dif- 
ferent peaks in Fig. g reflects the field-induced changes in 
the magnetic structure; a more extensive diffraction ex- 
periment is required to determine the high-field magnetic 
structure. To determine the high-field correlation length, 
the data have been fitted to an intrinsic 3D Lorentzian 
cross section convoluted with the experimental resolution 
function. For the correlation length along the a-axis, we 
obtain S,a = lOOOAzt 500A. The transverse resolution is 
too coarse to allow us to determine the correlation length 
in the transverse direction (0.5, 3k, k). 



The temperature dependences of the peak intensities, 
the separation between the incommensurate peaks, and 
the full peak width at H=6.9 Tesla are shown in Fig. |l^. 
The correlation length at this field is always finite, and 
no well defined phase transition is observed. Note that 
at high temperatures the period of the incommensurate 
spin modulation reverts back to the value characteristic 
of the zero-field magnetic structure. 

The unexpected result of our experiment is the finite 
correlation length of the high-field state. The domain 
size in this state (500-1000 A) is certainly large. How- 
ever, as the data of Fig. O show, there is no well de- 
fined phase transition as a function of temperature in 
the high-field state. This type of behavior is often found 
in magnetic systems with quenched disorder such as, for 
example, the doped spin-Peierls compound CuGeOs or 
Random Field Ising materials |^,Q . In some of these 
materials, the application of a magnetic field also results 
in the destruction of long-range order Ig^l- Therefore, 
one possible explanation of the finite correlation length 
in the magnetic field is that it is due to impurity effects. 
The magnetic long-range order can be destroyed, for ex- 
ample, by random fields that arise due to the presence of 
impurities in an antiferromagnetic system in which the 
applied field has a component along the direction of the 
spins ll^,^. In our samples, the loss of long-range order 
coincides with a transition from a state with magnetic 
periodicity which is a rational fraction of the unit-cell 
length (7/3) to a state which appears to be truly incom- 
mensurate. This gives an additional argument in favor of 
the important role of impurities because impurity effects 
are, in general, more pronounced in systems with con- 
tinuous symmetry than in systems exhibiting periodic 
"locked" structures. As mentioned above, the non-zero 
Curie constants found in our samples as well as our pow- 
der diffraction measurements suggest that either impuri- 
ties or structural disorder might indeed be present in our 
samples. 

While long-range order is found in the zero-field state, 
the very existence of magnetic ordering in CaCu203 may 
also be due to the presence of structural imperfections 
because the ground state of low-dimensional systems is 
extremely sensitive to even a small amount of disorder 
1^,^. It is well known, for example, that while pure 
SrCu203 stays in the spin- liquid state down to the low- 
est temperatures, substitution of just a few percent of 
copper ions with zinc induces low-temperature magnetic 
ordering 132]. At this stage, the precise control of the 
chemical composition of our samples has not yet been 
achieved. It is evident, however, that disorder can play 
a very important role in this material, and the synthesis 
of single crystals with a controlled amount of disorder is 
highly desirable. 



V. SUMMARY 

In summary, we have carried out magnetic suscepti- 
bility and neutron diffraction studies of the S=^ quasi- 
one-dimensional magnet CaCu203. Above T=50 K, the 
magnetic susceptibility of this material is in good agree- 
ment with the results of Monte Carlo calculations for 
a model consisting of ID spin-2 Heisenberg chains bro- 
ken into weakly-interacting segments of finite size with 
J||=2000±300 K, plus an additional small Curie term 
due to extrinsic impurities. CaCu203 undergoes a phase 
transition to an incommensurate magnetically ordered 
state at Tjv ~25 K. The ordering wavevector is (0.429(5), 
^, ^). The magnetic structure is incommensurate (or 
high oder commensurate) in the direction of the frus- 
trated interchain interaction. Weak commensurate mag- 
netic Bragg peaks with the reduced wavevector (i, i, i) 
are also observed below Tjv. The experimental data are 
consistent with a spiral (or possibly conical) magnetic 
structure with spins rotating in the ac crystallographic 
plane. The low-temperature magnitude of the ordered 
magnetic moment determined using this model magnetic 
structure is 0.2±0.07 fiB- 

Application of a magnetic field induces a metamag- 
netic transition at which the incommensurability of the 
magnetic structure is substantially reduced. Surprisingly, 
only a short-range order is found in the high-field state. 
The correlation length of the high-field state (500-1000 
A) is, however, relatively large. The finite correlation 
length found in this state presumably results from impu- 
rity effects, as discussed in section IV. 

In brief, CaCu203 is an experimental quasi-lD mag- 
netic system that provides an interesting opportunity to 
study the effects of quantum magnetism, magnetic frus- 
tration, and anisotropic magnetic interactions, and in 
which impurity effects can play a significant role. These 
subjects are currently attracting an intense attention of 
condensed matter physicists. Further experimental and 
theoretical investigation of CaCu203 is, therefore, of con- 
siderable interest. 

We would like to thank G. Shirane, A. Aharony, A. B. 
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TABLE L Atomic coordinates and temperature parame- 
ters for CaCu2 03 at T=10 K. The space group is Pmmn. 
The lattice parameters are a=9.9491(l) A, 6=4.0775(1) A, 
c=3. 4604(1) A. The atomic positions are given in the lat- 
tice parameter units. The temperature parameters B for the 
oxygen atoms were constrained to be equal in the refinement. 



Atom 



y 



B{k^) 



Ca 


0.25 


0.25 


0.3630(5) 


0.71(2) 


Cu 


0.0829(1) 


0.75 


0.8436(3) 


0.64(2) 


0(1) 


0.25 


0.75 


0.5832(4) 


0.94(2) 


0(2) 


0.0812(1) 


0.25 


0.8720(4) 


0.94(2) 



TABLE II. Observed integrated {lobs) and calculated in- 
tensities (Icaic) of incommensurate magnetic Bragg reflections 
in CaCu203. The measurements were performed at T=12 K. 



(h, k, 1) 


lobs 


Icalc 


(0.43, -0.5, -0.5) 


1770(30) 


1674 


(1.43, 0.5, 0.5) 


2030(70) 


4010 


(1.43, -0.5, -0.5) 


2330(30) 


4010 


(2.43, -0.5, -0.5) 


580(20) 


829 


(5.43, -0.5, -0.5) 


162(23) 


172 


(0.43, -1.5, -1.5) 


136(12) 


201 


(2.43, -1.5, -1.5) 


144(20) 


222 


(3.43, -1.5, -1.5) 


150(30) 


110 


(4.43, -1.5, -1.5) 


390(110) 


357 


(0.43, 0.5, 1.5) 


200(20) 


309 


(1.43, 0.5, 1.5) 


1290(60) 


1298 


(2.43, 0.5, 1.5) 


225(20) 


362 


(3.43, 0.5, 1.5 


130(23) 


168 


(4.43, 0.5, 1.5) 


676(80) 


544 


(0.43, 1.5, 0.5) 


746(25) 


481 


(-0.43, 1.5, 0.5) 


706(35) 


481 


(1.43, 1.5, 0.5) 


1290(175) 


1920 


(2.43, 1.5, 0.5) 


363(20) 


557 


(3.43, 1.5, 0.5) 


210(50) 


191 


(4.43, 1.5, 0.5) 


360(25) 


533 


(5.43, 1.5, 0.5) 


70(40) 


100 


(0.43, 2.5, 0.5) 


268(15) 


104 


(-0.43, 2.5, 0.5) 


270(40) 


104 


(1.43, 2.5, 0.5) 


525(15) 


414 


(2.43, 1.5, 0.5) 


170(60) 


HI 


(3.43, 1.5, 0.5) 


110(25) 


47 


(4.43, 1.5, 0.5) 


190(20) 
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FIG. 1. (a) The atomic structure of GaGu2 03. (b) The 
projection on the crystallographic a6-plane, and (c) on the 
ac-plane. Ga atoms are not shown. The dotted line indi- 
cates the structural unit cell. In (c) arrows correspond to the 
zero-field model magnetic spiral structure, and r) and x des- 
ignate the angles between the adjacent spins, as discussed in 
the text. 
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FIG. 2. Temperature dependence of the magnetic suscep- 
tibility of GaGu2 03 in a magnetic field of 1 kG oriented along 
the three major crystallographic axes, with the temperature 
independent part subtracted. The inset shows the inverse 
magnetic susceptibility after subtracting the T-independent 
part. The solid line is the result of a fit to the Gurie form. 




FIG. 3. (a) Magnetization versus magnetic field for H||a 
on ramping tlie field up and down, (b) Magnetization versus 
magnetic field along the three major crystallographic axes. 
The temperature is 5 K. The data for the magnetic field 
aligned with the b axis are shifted up by 4 emu/mol for clarity. 
The inset shows the temperature dependence of the transition 
field He for the magnetic field parallel to the a axis. 
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FIG. 4. High temperature magnetic susceptibility of 
CaCu2 03 measured in a magnetic field of 1 Tesla with the 
Curie part subtracted. The magnetic field points along the 
crystallographic b axis. The solid and dashed lines are the re- 
sults of Bethe ansatz calculations for the ID Heisenberg spin 
1/2 chain for several values of the intrachain magnetic cou- 
pling J|| . The sudden drop of the susceptibility above T=700 
K is due to sample decomposition. 
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FIG. 5. High temperature magnetic susceptibility of 
CaCu2 03 measured in a magnetic field of 1 Tesla aligned 
with the spin-chain (fe-)direction. The solid line is the re- 
sult of a Monte Carlo calculation for the finite-segment model 
described in the text with J||=1940 K. The inset shows 
the high-temperature magnetic susceptibility on an expanded 
scale. The dashed lines are the results of Monte Carlo calcu- 
lations with J||=1500 K and J||=2500 K. The sudden drop of 
the susceptibility above T=700 K is due to sample decompo- 
sition. 
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FIG. 6. Neutron diffraction scan along the (h, 0.5, 0.5) 
direction at T=12 K in sample 1. The solid line is the result 
of a fit to three resolution-limited Gaussian peaks as described 
in the text. 
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FIG. 7. The intensity of the (1.429, 0.5, 0.5) magnetic peak 
as a function of temperature in sample 1. The solid hne is the 
result of a fit to a power law as described in the text. The 
inset shows the temperature dependence of the intensity of 
the (1.5, 0.5, 0.5) commensurate magnetic peak. 
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FIG. 9. Magnetic field dependence of the peak intensity 
(a), the separation between the incommensurate peaks (b), 
and the full peak width (c) at T=10 K. The data were taken 
in the vicinity of the (0.5, 1.5, 0.5) position, as shown in Fig. 
8. The scanning direction was (h, 1.5, 0.5). The "left", "cen- 
tral" , and "right" keys in (a) refer to the lower, intermediate 
(commensurate), and higher h peaks in Fig. 8. 
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FIG. 8. Neutron diffraction scans along the (h, 1.5, 0.5) 
direction at T=10 K in various magnetic fields. Each scan is 
shifted along the y-axis by a constant offset value. The solid 
lines are results of the fits as described in the text. 



11 



H = 6.9 Tesia 



1000 p 
800 - 
600 - 
400 - 
200 - 



0.065 - 
-^0.060 - 
^0.055 - 
£ 0.050 - 
i 0.045 - 
0.040 - 
0.035 - 



left 



central 

n 

right 



0.15 
^0.14 
^0.13 - 

c 0.12 - 
I 0.11- 

o 0.10 - 
(^ 0.09 - 



b) 



c) 



a) 



s a 2 



• i 



10 



15 

T(K) 



20 



25 



FIG. 10. Temperature dependence of the peak intensity 
(a), the separation between the incommensurate peaks (b), 
and the full peak width (c) at H=6.9 Tesla. The scattering 
geometry and the keys in (a) are the same as in Fig. 9. 
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